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We present high-precision measurements of the Nusselt number TV as a function of the 
Rayleigh number R for cylindrical samples of water (Prandtl number a — 4.4) with a 
diameter D of 49.7 cm and heights L = 116.3, 74.6, and 50.6 cm, as well as for D = 24.8 
cm and L = 90.2 cm. For each aspect ratio V = D/L = 0.28, 0.43, 0.67, and 0.98 the data 
cover a range of a little over a decade of R. The maximum R ~ 10 12 and Nusselt number 
M ~ 600 were reached for T — 0.43 and D = 49.7. The data were corrected for the 
influence of the finite conductivity of the top and bottom plates on the heat transport in 
the fluid to obtain estimates of Moo for plates with infinite conductivity. The results for 
Moo and T > 0.43 are nearly independent of T. For T = 0.275 Moo falls about 2.5 % below 
the other data. For R < 10 11 , the effective exponent j e ff of Moo — N R lc ff is about 
0.321, larger than those of the Grossmann-Lohse model with its current parameters by 
about 0.01. For the largest Rayleigh numbers covered for T = 0.98, 0.67, and 0.43, j e ff 
saturates at the asymptotic value j = 1/3 of the Grossmann-Lohse model. The data do 
not reveal any crossover to a Kraichnan regime with 7 e // > 1/3. 



1. Introduction 

Understanding turbulent Rayleigh-Benard convection (RBC) in a fluid heated from 
below [Siggia(1994), Kadanoff(2001), Ahlers, Grossmann & Lohse (2002)] is one of the 
challenging and largely unsolved problems in nonlinear physics. An important aspect is 
the global heat transport that is usually expressed in terms of the Nusselt number 

M = QL/XAT . (1.1) 

Here Q is the heat-current density, L the sample height, AT the applied temperature 
difference, and A the thermal conductivity of the fluid in the absence of convection. 
A central prediction of various theoretical models [Siggia(1994), Grossmann & Lohse 
(2000), Grossmann & Lohse (2001), Grossmann & Lohse (2002), Grossmann & Lohse 
(2004)] is a relationship between M, the Rayleigh number R — ag AT X 3 / nv (a is the 
isobaric thermal expension coefficient, k the thermal diffusivity, v the kinematic viscosity, 
and g the acceleration of gravity), and the Prandtl number a — v/k. A model developed 
recently by Grossmann & Lohse (2000), based on the decomposition of the kinetic and 
the thermal dissipation into boundary-layer and bulk contributions, provided an excellent 
fit to experimental data of Xu, Bajaj & Ahlers(2000) and Ahlers & Xu (2001) for a 
cylindrical cell of aspect ratio T = D / L = 1 (D is the sample diameter) when it was 
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properly adapted [Grossmann & Lohse (2001), GL] to the relatively small Reynolds 
numbers of the measurements. However, the data were used to determine five adjustable 
parameters of the model. Thus more stringent tests using measurements for the same 
T but over wider ranges of R and a are desirable. A great success of the model was 
the excellent agreement with recent results by Xia, Lam & Zhou(2002) for much larger 
Prandtl numbers than those of Ahlers & Xu (2001), at Rayleigh numbers near 1.78 x 10 7 
and 1.78 x 10 9 . 

Here we present new measurements in a cell of diameter D = 49.7 cm for a ~ 4.4 that, 
for T = 0.98, extend to R ~ 10 11 . We also report results for D = 49.7 cm and T = 0.427 
and 0.667, as well as for D = 24.8 cm and T = 0.275. For T > 0.5 it is expected that 
the large-scale flow (LSF) in the cell [Krishnamurty & Howard(1981)] consists of a single 
convection roll [Verzicco & Camussi(2003)]. For T < 0.5, on the other hand, Verzicco & 
Camussi(2003) suggest that the system contains two or more rolls placed vertically one 
above the other. How this impacts the heat transport was one of the interesting questions 
to be addressed. Our results suggest that even the cell with T = 0.43 still contained only 
one roll because the data for Af fall on a smooth line drawn through those for the larger 
T. The T = 0.28 results fall about 2.5% below all the other data, suggesting a more 
complicated, perhaps two-roll, structure for the LSF. 

Most of our measurements were made at a mean temperature of 40°C, where a = 4.38. 
However, for T = 0.67 we also made measurements at mean temperatures of 50 and 30°C, 
corresponding to a — 3.62 and 5.42 respectively. We found a very gentle decrease of M 
with increasing a, approximately in proportion to cr~ - 044 . 

One of the experimental problems in the measurement of J\f(R) is that the side wall 
often carries a significant part of the heat current. Corrections for this effect arc not 
easily made, because of the thermal contact between the wall and the fluid that yields 
a two-dimensional temperature field in the wall, and because of the influence of lateral 
heat currrents through the wall on the fluid flow [Ahlers (2000), Roche et a/.(2001), Verz- 
icco(2002), Niemela & Sreenivasan (2003)]. The present project was designed to provide 
data that are not uncertain due to a significant side-wall correction. We used a classical 
fluid of relatively large conductivity confined by side walls of relatively low conductivity. 
The system of choice was water confined by Plexiglas with various heights L, and with 
the greatest diameter permitted by other constraints. We built four convection cells, 
three with D — 49.67 cm with heights L = 116.33, 74.42, and 50.61 cm, and one with 
D = 24.81 cm and L = 90.18 cm. For the V = 0.98 (L = 50.61) cell, which is most relevant 
to comparison with the theoretical model of GL, we estimated a wall correction [using 
Model 2 of Ahlers(2000)] of only 0.3% for R = 5 x 10 9 (A/" ~ 100) and smaller corrections 
for larger R. Based on this estimate we felt justified in neglecting the correction. 

A second experimental problem was pointed out recently by Chaumat et aZ.(2002) and 
by Verzicco(2004). Using direct numerical simulation, Verzicco(2004) showed that end 
plates of finite conductivity diminish the heat transport in the fluid when the Nusselt 
number becomes large. In a separate paper we shall give details of the apparatus used 
in our work [Nikolaenko et al. (2004)]. There we will describe measurements using two 
types of top and bottom plates of identical shape and size, but one set made of copper 
with a conductivity Ac« = 391 W/m K and the other of aluminum with A^; = 161 W/m 
K. That work yielded a correction factor that has been applied to the data reported here. 
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No 


T(°C) 


AT(°C) 


1Q- S R 


N 




No 


T(°C) 


AT(°C) 


10~ 8 i? 


TV 


Moo 


1 


39.932 


2.122 


576.8 


229.1 


229.5 


2 


39.964 


4.036 


1098.4 


281.4 


282.3 


3 


40.001 


5.935 


1617.2 


317.2 


318.6 


4 


40.044 


7.817 


2133.1 


346.6 


348.6 


5 


40.090 


9.689 


2648.5 


371.1 


373.6 


6 


40.141 


11.553 


3163.5 


392.4 


395.4 


7 


40.191 


13.405 


3677.1 


411.4 


415.0 


8 


39.997 


2.983 


812.7 


255.1 


255.7 


9 


40.032 


4.888 


1333.4 


298.5 


299.7 


10 


40.015 


3.934 


1072.6 


278.6 


279.4 


11 


40.012 


2.261 


616.4 


233.7 


234.1 


12 


39.990 


6.448 


1756.4 


325.5 


327.1 


13 


39.994 


13.406 


3652.1 


410.0 


413.6 


14 


39.973 


8.941 


2433.8 


361.2 


363.4 


15 


40.020 


10.809 


2947.2 


383.3 


386.1 


16 


39.998 


11.833 


3224.0 


393.7 


396.8 


17 


40.017 


15.710 


4283.1 


430.5 


434.7 


18 


40.007 


13.773 


3753.6 


413.2 


416.9 


19 


40.031 


17.636 


4810.5 


446.6 


451.3 


20 


40.002 


19.647 


5353.6 


462.0 


467.3 


21 


39.997 


9.872 


2689.6 


372.2 


374.7 


22 


39.997 


7.908 


2154.5 


347.4 


349.3 



Table 1. Results for r = 0.275 and D — 24.81 cm. In this and all following data tables two 
points are listed per line, and they are numbered in chronological sequence. 



No T(°C) 


AT(°C) 


HT 8 i? 


TV 




No T(°C) 


AT(°C) 


W- S R 


TV 


/Voo 


1 


39.849 


19.921 


11603.4 


595.4 


624.7 


2 39.976 


17.705 


10358.7 


574.5 


601.6 


3 


40.021 


15.661 


9177.2 


553.1 


577.9 


4 40.028 


13.693 


8025.7 


529.2 


551.7 


5 


39.993 


11.811 


6914.1 


505.3 


525.4 


6 39.996 


9.843 


5763.1 


476.9 


494.6 


7 


40.092 


7.689 


4517.2 


441.1 


455.8 


8 40.008 


3.936 


2305.5 


356.1 


364.8 


9 


40.009 


1.965 


1151.2 


286.0 


291.0 


10 40.053 


2.862 


1679.0 


322.4 


329.2 


11 


40.057 


4.824 


2830.3 


380.6 


390.8 


12 39.979 


5.958 


3486.2 


405.8 


417.8 


13 


40.014 


1.464 


857.9 


260.2 


264.1 
















Table 2. 


Results for Y — 


0.427 and D 


= 49.7 cm. 







2. Results 

2.1. The measurements 

Details of the apparatus and of experimental procedures are given by Nikolaenko et al. 
(2004). The precision of the measurements is typically near 0.1 percent, and systematic 
errors, primarily due to uncertainties of the diameter and length of the cell, are estimated 
to be well below one percent. Deviations from the Boussinesq approximation arc believed 
to be unimportant [Nikolaenko et al. (2004)]. We give the results in Tables 1 to 4. The 
measured TV derived from Eq. 1.1 as well as the Nusselt number Afx obtained after 
correction for the finite top- and bottom-plate conductivity [Verzicco(2004)] are listed. 
The relation between these is given by Af = /(X)TVoo where X is the ratio of the average 
thermal resistance of the plates to that of the fluid. We used the empirical function 
f(X) = 1 - exp[-(aX) b ], with the parameters a = 0.275 and b = 0.390 for D = 49.7, 
and a = 0.304 and b = 0.506 for D = 24.82, determined experimentally [Nikolaenko et al. 
(2004)]. The parameters were found to be independent of the aspect ratio but to depend 
on the plate diameter. 

2.2. Dependence on R and T 

The results for M x are shown on logarithmic scales in Fig. la and with greater resolution 
in the compensated form Moo/R 1 ^ in Fig. lb. We note that over most of the range of 
R the data for T > 0.43 reveal very little if any dependence of TVqo on T. Recent consid- 
erations by Grossmann & Lohsc (2003) had suggested a stronger T-dependence. Earlier 
experimental data [Wu & Libchaber(1992), Xu, Bajaj & Ahlers(2000)] also suggested a 
stronger T dependence; but those results were influenced by side- wall and/or end-plate 
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No 


T(°C) 


AT(°C) 


1Q- S R 


TV 




No 


T(°C) 


AT(°C) 


10~ 8 i? 


TV 




Moo 


1 


39 


.922 


17.639 


2693.9 


371 


.5 


389, 


3 


2 


39 


.816 


16.035 


2439.8 


359 


5 


376 


.0 


3 


39 


.923 


13.881 


2120.1 


343 


.9 


358. 


,8 


4 


39 


.937 


11.892 


1817.2 


327. 


7 


341. 


.0 


5 


39 


.984 


7.891 


1207.7 


288 





297. 


9 


6 


39 


.995 


3.955 


605.5 


231. 


5 


237. 


.3 


7 


40 


.013 


1.952 


299.1 


185 


.8 


189, 


1 


8 


39 


.979 


19.633 


3004.5 


384. 


2 


403 


.4 


9 


39 


.964 


20.637 


3156.4 


389. 


.9 


409. 


,7 


10 


39 


.898 


18.821 


2872.0 


378, 


4 


396 


.9 


11 


39 


.984 


16.693 


2555.0 


365 


.1 


382, 


2 


12 


39 


.915 


14.875 


2271.3 


351, 


4 


367 


.1 


13 


39 


.890 


12.974 


1979.2 


336 


.5 


350, 


6 


14 


40 


.109 


10.574 


1625.6 


316, 





328 


.3 


15 


39 


.974 


8.900 


1361.7 


299 


.4 


310, 


2 


16 


40 


.059 


2.846 


436.8 


209, 


.0 


213. 


.5 


17 


39 


.974 


17.693 


2707.1 


371. 


.7 


389, 


5 


18 


49 


.989 


19.566 


4101.7 


426, 





450. 


.1 


19 


49 


.886 


18.789 


3927.2 


420 


2 


443, 


6 


20 


49 


.969 


17.651 


3698.0 


412, 


,3 


434. 


.8 


21 


49 


.959 


5.003 


1047.9 


277 


.7 


286, 


.8 


22 


50 


.087 


6.710 


1410.6 


304, 


.8 


316. 


.0 


23 


50 


.136 


8.576 


1805.4 


329 


.9 


343, 


.4 


24 


50 


.146 


10.503 


2211.6 


351, 


2 


366 


.8 


25 


50 


.174 


12.406 


2614.4 


370. 


.4 


388, 


1 


26 


50 


.097 


14.525 


3054.1 


388, 


.4 


408. 


.1 


27 


49 


.902 


4.140 


865.7 


261 


.6 


269, 


.4 


28 


49 


.935 


3.092 


647.1 


239, 





245. 


.3 


29 


49 


.928 


2.120 


443.6 


212 


.4 


217, 





30 


49 


.973 


1.046 


219.1 


172, 


3 


175 


.0 


31 


39 


.974 


17.693 


2707.1 


371 


.7 


389, 


5 


32 


29 


.980 


19.647 


2007.0 


336 


1 


350 


.3 


33 


29 


.996 


17.656 


1805.0 


325 


.1 


338 


2 


34 


29 


.841 


16.006 


1624.5 


314, 


,4 


326 


.6 


35 


29 


.742 


14.246 


1439.0 


302 


2 


313, 


2 


36 


29 


.913 


11.963 


1218.2 


287. 





296 


.8 


37 


29 


.957 


9.918 


1012.1 


270 


.5 


278, 


9 


38 


29 


.932 


8.011 


816.6 


252. 


,8 


260. 


.0 


39 


29 


.859 


6.190 


628.8 


232 


.5 


238 


3 


40 


29 


.897 


4.157 


423.0 


205, 





209. 


.2 


41 


29 


.904 


2.171 


221.0 


167 





169 


5 


42 


30 


.036 


2.898 


296.8 


183, 


6 


186 


.8 


43 


30 


Oil 


4.913 


502.6 


216 


.6 


221, 


.5 


44 


30 


.018 


6.858 


701.8 


240, 


6 


246 


.9 


45 


29 


.992 


2.491 


254.6 


174 


.3 


177, 


1 





















Table 3. Results for T = 0.667 and D = 49.7 cm. 



No 


T(°C) 


AT(°C) 


10~ 8 i? 


TV 


Afoo 


1 


40.092 


1.792 


86.7 


125.3 


127.5 


3 


40.015 


3.898 


188.0 


158.8 


162.9 


5 


39.953 


5.969 


287.2 


181.2 


186.8 


7 


39.933 


7.952 


382.3 


197.9 


204.8 


9 


39.935 


11.829 


568.8 


224.2 


233.4 


11 


40.039 


15.500 


748.1 


244.3 


255.5 


13 


39.933 


17.645 


848.4 


254.1 


266.3 


15 


39.996 


1.489 


71.8 


118.0 


119.9 


17 


39.943 


17.633 


848.2 


254.2 


266.4 



Table 4. Results for T 



No 


T(°C) 


AT(°C) 


10~ f 


R 


TV 


TVoo 


2 


40.006 


2.940 


141. 


.7 


145.6 


148.9 


4 


39.975 


4.951 


238. 


.4 


170.9 


175.7 


6 


39.933 


6.979 


335. 


6 


190.2 


196.4 


8 


39.966 


9.832 


473. 


3 


211.4 


219.4 


10 


39.916 


13.803 


663. 


.3 


235.1 


245.4 


12 


39.407 


16.728 


789. 


6 


248.3 


259.9 


14 


40.038 


2.386 


115. 


1 


136.2 


138.9 


16 


39.935 


17.638 


848. 


2 


254.0 


266.3 



0.981 and D = 49.7 cm. 



effects. Note for instance that the side- wall correction made by Ahlers(2000) considerably 
reduced the V dependence originally seen by Xu, Bajaj & Ahlcrs(2000). Similarly, our 
present results for TV reveal some dependence on T, but the end-plate correction largely 
removes it. It is particularly noteworthy that the data for Y — 0.43 are not shifted signif- 
icantly relative to those for V = 0.67 because on the basis of the numerical calculations 
of Verzicco & Camussi(2003) one expects different structures for the LSF for these two 
cases. The data for T = 0.275 are lower by about 2.5%, suggesting that a transition in 
the LSF structure may occur between T = 0.43 and 0.275. 

A second important feature of the data is their dependence on R. Locally, over a limited 
range, the measurements can be fit by the powerlaw 



tV^ - N B?>'" 



(2.1) 
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Figure 1. (a): The Nusselt number Af x as a function of the Rayeigh number R on logarithmic 
scales, (b): The compensated Nusselt number Afoo/R 1 ^ 3 as a function of the Rayleigh number R. 
Open squares: F = 0.275. Solid circles: F = 0.43. Open circles: F = 0.67. Solid squares: F = 0.98. 
Solid lines: the prediction of Grossmann & Lohse (2001) for a — 4A. 



For the effective exponent we found j e ff = 0.323 for T = 1 near R = 2 x 10 10 and 
j e ff — 0.329 for T = 0.67 near R = 10 11 . As shown in Fig. 2, a single fit to most of 
the data for T = 0.98,0.67, and 0.43 yields j e ff — 0.322. All these values are close to, 
but definitely less than, the asymptotic large- R prediction 7 = 1/3 of the GL model. 
However, they are larger by about 0.01 or 0.015 than the GL prediction for T — 1 at the 
same R. This can also be seen qualitatively from Fig. lb where the GL prediction with 
its present parameter values is shown as a solid line. It remains to be seen whether the 
model parameters can be adjusted so as to reproduce this feature. 

Another interesting aspect noticeable in Fig. lb is that for T = 0.98, 0.67, and 0.43 
there is a seemingly sudden change in the i?-dependence of AAqo at large R to a power 
law with 7 = 1/3, i.e. with the asymptotic prediction of the GL model. This is illustrated 
more clearly in Fig. 2 which shows the relevant data with higher resolution. It is difficult 
to sec how the GL model can reproduce this rather sudden transition at finite R to 
its asymptotic exponent value. Rather, it seems likely that a new physical phenomenon 
not yet contained in the model will have to be invoked. The transition occurs for T = 
0.98, 0.67, and 0.43 near R = 5 x 10 10 , 9 x 10 10 , and 3 x 10 11 respectively and is reflected 
in the observation that the data for larger R fall on horizontal lines in the figures. In 
this range the measurements do reveal a dependence on T, with -/Vo(7e// — 1/3) = 
0.0606,0.0602, and 0.0595 for T = 0.98,0.67, and 0.43 respectively. One interpretation 
of these results for N is that the heat transport is diminished, albeit only very slightly, 
by a larger travel distance, and presumably a larger period, of the LSF. At smaller R 
an effective powerlaw with N = 0.0797, j e ff = 0.322 (solid line in Fig. 2) fits the data 
at all three T values within their statistical uncertainty, showing that the results are 
essentially T independent. It is a surprise that the data for V = 0.67 and 0.43 do not 
differ more from each other because different structures for the LSF had been predicted 
for these two cases. [Verzicco & Camussi(2003)] The T = 0.275 results do indeed have 
a Nusselt number that is smaller by about 2 to 3%, suggesting that the transition in 
the flow structure occurred between T = 0.43 and 0.275. Those result do not show the 
saturation at large R with -f e ff — 1/3, and lead to j e ff — 0.321 over their entire range. 
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CO 
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Figure 2. The compensated Nusselt numbers A r 00 /-R 1 ^ 3 as a function of the Rayleigh number R 
for F = 0.98 (solid squares), F — 0.67 (open circles), and T = 0.43 (solid circles) on an expanded 
scale. Lines: powerlaws A/"oo = NqR 1 " 11 . Solid line: iVo = 0.0797, y e ff = 0.3222. Dashed lines: 
7e// = 1/3 and (from top to bottom) 7V = 0.06063, 0.06019, 0.05945. 
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Figure 3. (a): The compensated Nusselt numbers Afoc/R 1 '' 3 as a function of the Rayleigh 
number R for T = 0.667 at three values of the Prandtl number, (b): The reduced Nusselt 
numbers Afaocr ' 044 / R 1 ^ 3 as a function of 7?. Up pointing triangles: 50°C and a = 3.62. Open 
circles: 40°C and a — 4.38. Down-pointing triangles: 30°C and a — 5.42. 



2.3. Dependence on the Prandtl number a 

The Nusselt number JVoc(a) has a broad maximum near a ~ 4. Thus the dependence of 
Af on a is very weak and difficult to determine from measurements with various fluids 
of different a because of systematic errors due the uncertainties of the fluid proper- 
ties [Ahlers & Xu (2001), Xia, Lam & Zhou(2002)]. We determined Mia, R) with high 
precision over a narrow range of a by using the same convection cell and by chang- 
ing the mean temperature. In that case errors from different cell geometries are largely 
absent and the properties are known very well. Similar measurements over the ranges 
5 x 10 6 < R < 5 x 10 8 and 4 < a < 6.5 were made by Liu & Ecke(1997). Their results 
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Figure 4. The reduced Nusselt numbers Afoccr ' 044 /R 1 ^ 3 as a function of R. (a): Stars: Data of 
Roche et a£.(2004) for 2.5 < a < 6.0 and T = 0.50. Other symbols: this work for a = 4.4 and 
r = 0.98 (solid squares), 0.67 (open circles), 0.43 (solid circles), and 0.28 (open squares), (b): 
open diamonds: data of Chavanne et a/. (2001) for T = 0.5 and 2.0 < a < 4.5; solid symbols and 
open circles: this work for T = 0.98, 0.67, and 0.43. 



can be represented well by 

N = N oa cr a ^" K 7 '" (2.2) 

with 7 e // = 0.286, a e // = —0.030, and iV o = 0.1780. The negative value of a e ff 
indicates that the maximum of N(a) is below a ~ 4. 

We used our T = 0.67 cell and made measurements at 30, 40, and 50°C corresponding 
to a = 5.42, 4.38, and 3.62 respectively. The results are included in Table 3, and in the 
compensated form jV^/i? 1 / 3 they are shown as a function of R in Fig. 3a. Over the range 
2 x 10 10 < R < 4 x 10 11 they yield a e // — —0.044, indicating that also for this range of R 
the maximum of M{a) occurs below a = 4. The reduced Nusselt number A/" oo o' 044 /i? 1/ ' 3 
is shown in Fig. 3b. It shows that the data, within their experimental uncertainty, collapse 
onto a unique curve when divided by o~ a ° s f . The observed cr-dependence is somewhat 
stronger than that of the GL model with its present parameter values. 

2.4. Comparison with other results 

Measurements oiJV(R) for T = 0.50 over the range 10 8 < R < 5 x 10 10 and a wide range 
of a were made at cryogenic temperatures using gaseous helium by Roche et al.(2QQA). 
A direct, highly quantitative comparison with our results is possible only for the data 
points with cr-values fairly close to ours where the cr-dependence of N can reasonably 
be expected to be given by Eq. 2.2 with a e ff ~ —0.044. In Fig. 4a we show results of 
Roche et al. (2004) for A r oo o' 044 /i? 1 / 3 as a function of R. These data were corrected by 
the authors for the side- wall conductance, using a procedure described by them. Because 
of the small conductivity of helium gas one expects end-plate corrections to be negligible 
in this case. We also show our results for comparison. One sees that the data of Roche 
et al. (2004) fall into three well defined groups, with a nearly uniform vertical spacing 
between them of close to four percent. There is excellent agreement /consistency of the 
uppermost branch with our data for T = 0.98, 0.67, and 0.43. This is consistent with 
the absence of any significant aspect-ratio dependence in our data. The existence of the 
lower two branches is more difficult to reconcile with our results. Roche et aZ.(2004) 



8 Alexei Nikolaenko, Eric Brown, Denis Funf schilling, and Guenter Ahlers 
[2 1 1 1 1 A 1 1 I 



A 




R 

Figure 5. Compensated Nusselt numbers Moo/ 7? 0,309 as a function of 7?. Open circles: Niemela & 
Sreenivasan (2003), F = 1, 0.7 < a < 0.9. Open triangles: Niemela et al. (2000), F = 1/2, a ~ 0.7. 
Open diamonds: Chavanne a/.(2001), T = 1/2,0.7 < a < 2.0. Solid circles and stars: present 
work, 0.43 < V < 0.98, a = 4.4. Solid squares: present work, N X ,F = 0.28, a = 4.4. Dash-dotted 
line: powerlaw fit to our data with F > 0.43 (7e// = 0.3207). Dotted line: powerlaw fit to our 
data for F = 0.28 {~/ eff = 0.3193). Solid (dashed) line: GL model for a = 4.4 (0.8). 

suggest that they encountered more than one distinct state of their LSF. In our work 
we never found multi-stability for any T, and the data for T = 0.67 and 0.43 (which 
according to the calculations of Verzicco & Camussi(2003) should correspond to different 
flow structures) agree with each other and are consistent with the upper branch of Roche 
et al., at least in the range where Jeff has not yet saturated at "feff — 1/3. Our data for 
T = 0.28 are lower than those for our larger T, suggesting a difference in the LSF and 
indicating that the transition from a single cell to a more complicated structure occurs at 
r < 0.43 in our system. However, our results for T = 0.28 are only about 2 or 3 % lower 
than the larger-r data and not as low as the results from the middle or lower branch of 
Roche et al. 

In Fig. 4b we compare our results with those of Chavanne et aZ.(2001) from cryogenic 
experiments for 2.0 < a < 4.5 and T = 0.50. The end-plate corrections are expected to 
be negligible. Because of the low conductivity of the fluid, the side-wall contribution to 
the conductance of the cell is significant, but apparently no correction was made. The 
authors interpret their data to imply j e ff ~ 0.38, and attribute the large value to a 
breakdown of the boundary layers adjacent to the top and bottom plates. In that case 
one expects that a new regime first proposed by Kraichnan(1962) with an asymptotic 
exponent 7=1/2 (and logarithmic corrections) should be entered. Our data do not 
reveal such a large exponent, and in the overlapping range 10 11 < R < 2 x 10 12 remain 
consistent with J e ff = 1/3. 

Finally, in Fig. 5 we compare our results for T = 0.98, 0.67, and 0.43 (solid circles 
where j e ff ~ 0.32 and stars where j e ff ~ 0.333) and T — 0.28 (solid squares) on a less 
sensitive vertical scale with several measurements at cryogenic temperatures. A more 
comprehensive comparison was presented by Niemela & Sreenivasan (2003) (NS) [we 
choose the same representation in terms of AA00/-R ' 309 that was used by them in their 
Fig. 5]. There are systematic differences between the data sets that are, according to 
NS, larger than possible systematic errors in the experiments. Our data show very little 
change of -f e ff with R, and the solid circles in Fig. 5 yield 7 e // ~ 0.321. The data by 
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NS (connected open circles) show a dependence on R that differs from ours, with 7 e // 
varying from about 0.28 to about 0.35 as R changes from 10 9 to 10 12 . The data of Niemela 
et al. (2000) do not show such a strong variation of 7 e // and can be described well by 
a single 7 e // = 0.30 over the i?-range of the figure. The data of Chavanne et oZ.(2001) 
(open diamonds) are on average slightly larger than ours but have an independence that 
is consistent with that of ours. 

We have been unable to rationalize the large variations of the effective exponents of 
Af from one experiment to another, and in the case of the data of NS within a single 
experiment with R, in terms of known fluid mechanics. NS suggested that differences in 
the LSF structure are responsible, and state that "such differences seem to arise from 
delicate interplay among detailed geometry as well as Prandtl and Raylcigh numbers". 
Looking at our data, the absence of any significant T-dependence is apparent from the 
clustering of the data for T > 0.43 (solid circles) near the straight dash-dotted line 
corresponding to j e ff — 0.3207; i.e. T per se does not play a major role. Even for 
T = 0.28, where the data are displaced vertically in the figure, presumably because of a 
change in the LSF, they have nearly the same effective exponent 7 e // = 0.3193 (dotted 
line) . Thus it appears from our data that differences in the LSF do not have a large effect 
on the independence of A/^i?). We wish we had a more satisfying explanation of the 
different results obtained by the various experiments. 
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